Quantum master equations from classical Lagrangians with two stochastic forces 
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We show how a large family of master equations, describing quantum Brownian motion of a har- 
monic oscillator with translationally invariant damping, can be derived within a phenomenological 
approach, based on the assumption that an environment can be simulated by two classical stochas- 
tic forces. This family is determined by three time-dependent correlation functions (besides the 
frequency and damping coefficients), and it includes as special cases the known master equations, 
whose dissipative part is bilinear with respect to the operators of coordinate and momentum. 
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I. INTRODUCTION 

The phenomena of irreversibility and damping in quan- 
tum systems have been the subjects of numerous studies 
from the very first years of Quantum Mechanics (QM) 
P, 0, H, HI- They are attributed to the action of some 
environment that drains information from quantum sys- 
tem iS, thus increasing the entropy during its evolu- 
tion. The environment is usually referred as reservoir 
1Z whose variables are expressed as g-numbers (opera- 
tors acting in the Hilbert space), and it is assumed to 
have an infinity of degrees of freedom. In the most gen- 
eral case, one has to treat the whole quantum system 
S + 1Z, taking into account all details of each subsystem 
and the interaction between them. It is known, however, 
that under certain conditions the influence of environ- 
ment can be described approximately by means of a few 
parameters which enter some dynamical equations con- 
taining only the variables related to the system S. In 
the Schrodinger picture such equations for the statisti- 
cal operator p{t) or its different representations (density 
matrix, Wigner function, etc.), known under the name 
quantum master equations (^QMEJl, were studied in nu- 

"01110 
I22L I22J, 123|, |24|, |25(. An alter- 
native description in the Heisenberg picture is achieved 
within the frameworks of quantum Langevin equations or 
Langeyin-Heisenberg equations (LHE's) flri 0, [13, [H, 
HI Mi, SH HI |M IS M HIES, GUT These two 
approaches were discussed in detail in many books and 

reviews ESEiiiiiiiisiiiiiiiiiiiiiialiSIHilil. 

The shortest and simplest formal way of introducing 
relaxation in quantum mechanics is to postulate some 
general structure of the master or Langevin-Heisenberg 
equations, which guarantees fulfilment of the fundamen- 
tal requirements of QM (such as the properties of her- 
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miticity and positivity of the statistical operator, as well 
as its normalization; or preservation of canonical commu- 
tation relations between the time-dependent operators) 
for any instant of time. For example, a general struc- 
ture of dynamical mappings of density operators preserv- 
ing positive semidefiniteness was established in H, H3] , 
and equivalent differential equations were considered in 
[Hfil Ha, H3, HI] . For a more restricted problem of quan- 
tum Brownian motion this approach was used, e.g., in 
Refs. [H, [6(J EH, E3, [H, [64|, where some sets of free 
parameters were chosen in such a way that mean val- 
ues of coordinates and momenta obeyed given classical 
equations of motion. 

However, it seems desirable to have some schemes 
permitting to derive master equations from some gen- 
eral principles. The most natural way to do this is to 
start from some explicit Hamiltonian for the total system 
S + 1Z. Then the dynamical equations for the subsystem 
S arise as a result of taking the trace over the reser- 
voir variables. In this approach both the subsystem and 
reservoir are considered as quantum objects from the very 
beginning. A frequently used simple model of the reser- 
voir consists of an infinite number of harmonic oscillators 
(HO's) with a frequency distribution chosen according to 
some hypothesis. Note, however, that although the use of 
a reservoir is usually viewed as a "microscopic" approach 
for the derivation of LHE or QME, in reality, unless one 
knows precisely the nature of the environment and its 
interaction with S (described by a Hamiltonian with pa- 
rameters dependent on fundamental constants; see, e.g., 
[65[ for this exceptional case), strictly speaking, it should 
be considered as phenomenological. This occurs because 
several assumptions must be made concerning the very 
nature of the reservoir modes, on the parameters of the 
reservoir Hamiltonian and on the 1Z — S interaction term. 

Taking this point of view, it seems interesting to study, 
what kinds of master equations can be obtained, if one 
follows the phenomenological path from the very begin- 
ning to the very end. This is the aim of our paper. As a 
general principle we assume the Lagrange-Hamiltonian 
formalism, which impelled the development of physics in 
general and in particular QM and field theories. It is 
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complemented by the assumption that S interacts with 
time-dependent (TD) classical fluctuating forces [-Fj(i)], 
which simulate the action of the environment without 
any further preoccupation about the nature of their mi- 
croscopic origin. We show that combining some phe- 
nomenological Lagrangian, which takes into account the 
effect of dissipation, with a suitably chosen set of classi- 
cal stochastic forces and following standard rules of the 
"canonical quantization" procedure, one can obtain in a 
quite straightforward and simple manner a large class of 
master equations, which embraces many equations, intro- 
duced earlier within the frameworks of different schemes. 
This approach was used earlier in Refs. [H, H3, US IsS [zil • 
However, in those papers the authors considered only one 
classical stochastic force. As a result, the equations ob- 
tained had some restrictions on their domain of validity. 
We show that introducing two stochastic forces one can 
obtain more general equations. 

The paper is organized as follows. In Sec. [IT] we give a 
brief review of known master equations describing quan- 
tum damped oscillators. New results are contained in 
Sec. IIII1 where we derive the master equation with the 
use of two effective classical fluctuating forces. We show 
that in the classical case, it is always possible to use 
a single force, without changing the physics. However, 
in the quantum case two forces are necessary to obtain 
equations compatible with the principles of QM. Sec. IIVI 
contains a summary and conclusions. In appendix A we 
expose some details of calculations. 



II. CONVENTIONAL MASTER EQUATIONS 

The Hamiltonian Hq considered in this paper corre- 
sponds to a particle of mass m and angular frequency luq 
subjected to a one-dimensional harmonic force, 

H a = Hluq (a f a + 1/2) = p 2 /(2m) + mt^i 2 /2, (1) 

where the non-hermitian lowering operator a can be writ- 
ten in terms of hermitian operators, the particle position 
x and its momentum p, as 

a = (2ftwo m ) (mujox + ip) . 

The general form of the master equation is 

dp /alt + (i/h) [H, p] = Cp, (2) 

where the dissipative nonunitary superoperator C ac- 
counts for the influence of the environment and the op- 
erator H can be different from Hq. One of the most fre- 
quently used master equations for the damped harmonic 
oscillator has an origin in the problems of quantum op- 
tics. It corresponds to H = Hq and 

Cp = X (n + 1) (2dpa ! — a) dp — para) 

+Xn (2a) pd — da) p — pad)) , (3) 



where fi = n(ujQ,T) — (e 13 — 1) _1 is the reservoir mean 
number of quanta at temperature T (j3 — /iwo/fcs? 1 ), 
ks is Boltzmann constant and A is the decay rate or 
damping constant, which contains the strength parame- 
ter of the interaction between 1Z and S. The operator 
§5§ seems to be derived for the first time under the as- 
sumption of weak coupling between a quantum system S 
and a quantum thermal reservoir and by adopting the 
Born-Markov approximation in [IH, [H, [TJ [H| . Its rep- 
resentation in the Fock basis was derived in the general 
form in d, [t| [nj EH, although the special case of the 
equation for diagonal matrix elements can be traced to 
papers by Landau, Bloch and Pauli [H,[l,[l] (see also [(|). 
An immediate generalization of operator is 

Cp = £ (2$^ - - P&&) , (4) 

3 

where $j are arbitrary linear operators (their number 
also may be arbitrary). Sometimes the right-hand side 
of is called the "Lindblad form" , after the study [58| , 
although this general structure, preserving the hermitic- 
ity, normalization and positivity of p, was discovered by 
several authors earlier [55|, [56|, [13] • 

The relaxation operator in terms of operators x and p is 
usually associated with the problem of quantum Brown- 
ian motion. The most general master equation, preserv- 
ing the normalization and hermiticity of the statistical 
operator p and containing only bilinear forms of opera- 
tors x and p, corresponds to the choice 

H = Hq + ^{x, P }, (5) 



Cp = 2h ^' ^' ^ ~ 2h ^' ^' ^ ~ ~W ^' ^' ^ 
~%~ fa fa ft] + Tpr & fa + "jpr fa t f ' ^ ^ 

where p, A, D x , D p and D z can be, in principle, arbi- 
trary functions of time. This form was considered, e.g., 
in Refs. [3(| l46l l6lll (in the case of time- independent 
coefficients) and [3 71 ] (with arbitrary time-dependent co- 
efficients). Some differences in explicit expressions can 
be removed with the aid of identities [x,p] — ih and 

[x, {p, P}} + \p, {x,p}} + [p, {p, x}] = 0. 

The meaning of parameters A and p becomes clear from 
the equations for mean values of the coordinate and mo- 
mentum: 

d(x)/dt = (p)/m+ (p - X)(x), (7) 
d(p)/dt = -mu%(£) - (p + X)(p). (8) 

The choice A = p eliminates the friction term from the 
equation for d(x)/dt. This special case was studied by 
Dekker in [68| . Generalizations of ([5]) and ^ to multidi- 
mensional systems (in particular, to the case of a charged 
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particle in a magnetic field) were given in [59|, [6(j ■ The 
superoperator §5§ is a particular case of ^ for 



Y, a *iPi =iX/(2h)-D z /h 2 



(17) 



D 



X 



p - -mhwo I n + - 



1 



(mLOo) 2 D x 



fi = 0. 



The operator master equation being written in 
terms of the Wigner function, 

W (x,p, t) = ^J dye- 2 ^/ h (x ~ y\p(t)\x + y) , (9) 

assumes a simple form of the Fokker-Planck equation, 



^-(\p/m + (fi-X)x] W) 



d 2 W 
dp 2 



D 



d 2 W 
~dx T 



2D. 



d 2 W 
' dpdx 



(10) 



Thus we see that the terms proportional to D x , D p , and 
D z in Eq. ^ describe the diffusion in the phase space. 

Introducing the phase space vector variable q = (x,p), 
one can write Eq. (|10p in a compact form 



dW (q, t) 
dt 



d d 2 W 

Aq^+Aiff- 

% dq^dqj 



(sum over repeated indices is understood), where 



/i — A 



—mtu 2 , —(fi + X) 



(11) 



(12) 



is a 'drift' matrix, which governs the evolution of the first 
order statistical moments (mean values), 



d(q)/dt = A(q). 



(13) 



Introducing the covariances djk = \{qj<j_k + — 
(<7j) (<7fc), one can verify that both equations, (|10p and 
(jlip , result in the following equation for the symmetrical 
covariance matrix M = 1 1 cr ^ ^ 1 1 : 



dM/dt = AM + MA + 2D, 



(14) 



where A is the transposed matrix of A and D = ||-Dy|| 
is the symmetrical diffusion matrix (D12 = D 2 i = D z ). 

Taking the operators <tj in Eq. (UJ) as linear combina- 
tions of operators x and p, 



<&j = ajX + f3jp 7 



(15) 



one can verify that the operator ([6]) can be rewritten in 
the form (H|), provided the following conditions can be 
satisfied [71(: 



5>/ = zvn 2 , y,\ 



D p /h 2 



(16) 



In view of the Schwartz inequality, 



Ei^-i 2 Ei 



> 



the condition of compatibility of Eqs. (JT6J) and (jTTJ) is 
the inequality 



D P D X -D 2 Z > (HX/2) 2 = [ftTr(A)/4] 



(18) 



The condition (|18p , which was derived and discussed from 
different points of view in Refs. [13, M, HI HI, M, M, 
[721 . [T3l ] , guarantees that the positivity of statistical oper- 
ator is preserved for all times and for any physically ad- 
missible initial state. This is the necessary and sufficient 
condition (together with conditions D x > and D p > 0) 
of reducibility of the operator ((6|) to the Lindblad form 
dU [59j, [60j . Note that parameter p does not enter the 
constraint (fT8|) . because it is related to the correction to 
the Hamiltonian §5§ and not to the non-unitary part of 
the total Liouville superoperator. So its presence in the 
"friction forces" in Eqs. ([7]) and ([5]) is not relevant to the 
existence or non-existence of the Lindblad representation 
of the master equation. 

Some frequently considered master equations with 
time-independent coefficients, such as, e. g., the Agar- 
wal equation [74| (with X = fjt, D p = ImXiOtfh and 
D x = D z = 0) or its special case, known as the Caldeira- 
Leggett equation [751 ] (with A = /i, D p = 2mAfcsT 
and D x = D z = 0), do not satisfy the condition (|18p . 
Consequently, these equations can result in violations 
of the positivity of the statistical operator (which is 
equivalent to the violation of the uncertainty relations 
[47I I60L l70l . [7l| ) at the intermediate stages of evolution, 
if one tries to use them outside the domain of their valid- 
ity (which corresponds to the limit of high temperatures, 
k B T » huj ). 

On the other hand, it was shown in [zl, [z3, [zl, [zl, [H 
that the Wigner function of a subsystem S interacting 
with a reservoir 1Z satisfies Eq. (jTTJ) at any time, if, (I) 
the total Hamiltonian of the whole system S + 1Z is an 
arbitrary quadratic form with respect to coordinates and 
momenta (in particular, the interaction Hamiltonian can 
be an arbitrary bilinear form with respect to the coordi- 
nates of S and 1Z) and if, (II) the initial statistical oper- 
ator of the total system is factorized, p tot = p s p-ji, where 
Pn is an arbitrary Gaussian state (i.e., not necessarily 
thermal, it can be squeezed, for example). However, in 
such a case, (a) the matrices A and D are explicitly time- 
dependent (which is interpreted sometimes as a manifes- 
tation of non-Markovian evolution [371. l8ll \82\), and (b) 
the constraint (|18p (or its multidimensional generaliza- 
tions [l^HO]) can be violated. This does not mean that 
the state Ps(t) can become unphysical - no, simply in this 
case the system S does not pass over all possible mixed 
states in the process of evolution, but it moves only along 
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some specific trajectories in the Hilbert space; see in this 
connection also [36|, IzE HH • 

Time-independent diffusion and drift matrices appear 
only asymptotically, as t — > oo (physically, after some 
characteristic time determined by the properties of TV). 
Moreover, the sets of diffusion coefficients satisfying (fl8|) 
can be obtained only for specific forms of the 1Z — S in- 
teraction Hamiltonian. For example, in the case of a 
thermal reservoir, the relaxation superoperator ([2]) can 
be derived if the interaction only has the so called rotat- 
ing wave approximation form, which results in the drift 
matrix (TT2"]) with [i = [79j]. The superoperator ([5]), with 
arbitrary time independent diffusion coefficients, can be 
derived with the use of squeezed (or rigged) reservoirs 
[33L HH . Our main goal in this paper is to find the sets of 
drift and diffusion coefficients (possibly, time-dependent) 
that can be obtained from the scheme of quantization of 
classical equations with two stochastic forces. 



III. FROM CLASSICAL STOCHASTIC FORCES 
TO QUANTUM MASTER EQUATIONS 

A. Classical treatment 

A typical equation of motion of a particle of mass to 
subjected to a linear friction force in one dimension is 



x + r t i + 



1 dV(x,t) 
m dx 



0. 



(19) 



where V(x, t) is a potential and T t a TD dissipative func- 
tion. Although Eq. (fill)]) describes a nonconservative sys- 
tem, it is known for a long time (see, e.g., [H, H|| and 
references therein) that it can be derived from some La- 
grangian. The most simple one is known under the name 
Bateman-Caldirola-Kanai Lagrangian [13, , 



L (x, x, t) 



1 



-mx 



V(x,t) 



(20) 



We consider a simple generalization of (|2U)l . 
1 



L (x, x, t) 



-mx 2 -V(x,t) + xF t + xG t )e u , (21) 



where Ft and Gt are arbitrary generalized TD forces as- 
sociated to the position x and to the velocity x. The 
canonical momentum is 



P = dL/dx = e r * {mx + G t ) 



(22) 



and we define the physical momentum as 

p = mx + G t =Pe- Tt . (23) 
The Hamiltonian associated to the Lagrangian (f2"Tj) is 



p2 

Hit) = Px-L = —e- Tt + [V(x, t) - xF t ] e Ft 
2m 

-PGt/m + G\ exp(r t )/2m. (24) 



The classical equations of motion for canonical coordi- 
nates (Hamilton equations) are 

x = dH (t) I BP = (Pe- Tt - G t ) /to, 

P = -dH(t)/dx = {-dV(x,t)/dx + Ft)e Tt , 

whereas for the physical coordinates one obtains 

x = p/m — Gt/m, p — —dV(x, t) /dx + F t — pT t , 

or as a single second order equation (Newton equation) 

1 dV{x,t) 1 



x + T t x 



m dx 



-(Ft- T t G t -Gt), (25) 



where the RHS contains TD terms only. 

We see that in classical mechanics, where the coordi- 
nate x is the only independent variable (since p or P 
are functions of x), the presence of two terms, —xF t e Tt 
and —PGt/m, in Hamiltonian (|24|) is redundant, be- 
cause the dynamics depends only on the combination 
T{t) = F t -t t G t - G t . A usual choice is G t = 
and F t = T{i). But one can obtain the same dynam- 
ics, choosing F t = and finding the function Gt from 
the equation — (T t G t + G t ) — F{t), whose solution is 



G t = Ke~ Tt 



,-r, 



where K is an arbitrary constant. 

However, both the forces, F t and Gt, are important 
in the quantum case, because of the non-commutation 
property of position and momentum. These forces give 
different contributions to the dynamical evolution of the 
system state, as we show in the following subsection. 



B. Quantum treatment 



Having the Hamiltonian function (|24p . one can try to 
"quantize" the classical dissipative system, transforming 
(I24p to an operator by means of the usual rules and writ- 
ing the time-dependent Schrodinger equation with this 
"quantum" Hamilton operator. This idea was formulated 
for the first time by Caldirola and Kanai (CK) in 1940s 
[88j , and since that time it was developed or criticized by 
many authors (see, e.g., Refs. [4(|[89j]. It is known by 
now that such a simplified approach suffers from many 
drawbacks. For instance, the CK Hamiltonian is explic- 
itly time-dependent, so it is closer to the system with 
time-dependent mass than to the genuine dissipative sys- 
tem. Moreover, the problem of finding the Hamiltonian 
for the given equations of motion has no unique solution, 
and pra ctically all such Hamiltonians have some pathol- 
ogy [H, But the main physical defect of the CK 
scheme is that it implies that the quantum state of the 
system remains pure during the evolution, because reg- 
ular classical fields interacting with a quantum system 
S do not change its informational content as time goes 
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on, even if energy is not conserved; an initial pure state 
■0,5(0) will evolves as a pure state ^Jg(t). On the other 
hand, it is known that dissipation is connected with a loss 
of quantum purity. Thus one has to describe the system 
in terms of the density matrix or its equivalent forms, 
such as the Wigner function, for example. But how to 
find equations of motion for the density matrix? 

An answer was given in Refs. [66|, [6?], [H, [6^, l70l |: 
one should start not from the Schrodinger equation for 
a wave function, but from the von Neumann-Liouville 
equation for the statistical operator, considering F t and 
Gt as stochastic forces and performing some averaging 
over these forces. This averaging results in an informa- 
tion drain from the quantum system S. So irreversibility 
is verified and the entropy of the system changes due to 
the random character of the classical fields. However, 
by using a single stochastic force one obtains equations 
which not always preserve the property of positivity of 
the statistical operator. Our goal here is to show that 
correct equations, satisfying all principles of QM, can be 
derived in a very simple way, if one introduces two classi- 
cal stochastic forces, instead of a single one, besides one 
regular dissipative function. Such a derivation was not 
done earlier, as far as we know. 

The equation for the time evolution of the density ma- 
trix following from Hamiltonian (|2"4")l reads 



dt 



1 



H ef (t) + W(x,P,t),p t 



where 



Hef (t) 



t P 2 /(2m) + V(x,t)e 



(26) 



(27) 



Iterating Eq. (|32[) and deriving with respect to time we 
get the equation 



d Pt 
dt 



' [ W(t), Po 



ih 



(ih) Jo 



dt' 



W(t), W(t'),p t > (33) 



which is still exact. Its physical content is the same as 
Eq. (J30J) . Using the RHS of (32]) for p t , and inserting it 
in (|33p recursively, we obtain an infinite series 



dp t 
dt 



1 



1 



k=l 



W(t), W(t 



(ih) 

[w(t 



dtx 



k),P 



dth 



(34) 



The next step is to perform averaging over stochastic 
forces in Eq. (|34|) . We assume that classical stochastic 
functions Ft and Gt are Gaussian with zero average val- 
ues, F t = Gt = 0, and that they are delta-correlated with 
a TD function (not a stationary process): 



F tl F t2 = 2A tl 6(t 1 -i 2 ), 
G^G7 2 = 2m 2 B tl 6{t 1 -t 2 ), 
F)~G7 2 = 2mC t J(h-t 2 ). 



(35) 
(36) 
(37) 



Moreover, we take into account the important 
properties of Gaussian stochastic processes, namely, 
J tl J* 2 ... Jt 2n+ i = for an odd number of terms (J t stands 
for F t or Gt), whereas for an even number 



W(x,P, t) 



l xF t - PGt/m. 



(28) 



W(x,P, t) is a stochastic operator if F t and Gt are as- 
sumed as stochastic forces. The TD term G\ exp(r t )/2m 
was thrown out from H e f (t) since it does not contribute 
to the equations of motion. 

Using the unitary evolution operator XJ t corresponding 
to the effective free Hamiltonian (|2"T1) . 



ihdU t /dt = H ef (t)U t , 
we make a unitary transformation 

Pt = u t p t u}, 

which removes the term H e f (t) from Eq. 



ihd~p t /dt= \ W{x t ,Pt,t),~Pt\ = [W(t),~Pt 
and where 

x t = UjxUt, P t = u\PU t . 
A formal solution to Eq. ([ST))) is 



Pt = Po + 



4 I dl ' 

in 



W(t'),~p t , 



(29) 
(30) 
(31) 
(32) 



J tl Jt 2 ---Jt 2n — ^ JtiJtj ■ Jt k Jt n (38) 

all pairs 

where the average is over ensembles. Since stochastic 
operators W(tk) are linear combinations of i tfc and Pt k , 
only the terms with even numbers of operators W(tk) 
survive after the averaging in Eq. (|34[) . so that we arrive 
at the series containing only even powers of H: 



dp t _ \ - ' ! - - M 

dt ~ Aj(ift)2k*»w 



(39) 



The first term of this expansion is a sum of four integrals 
containing double commutators, 



dt'e Tt+T <'F t F v [x u [x v ,p Q \] 



+m 1 e Tt 



dt'FtGv 



dt'e r t'G t F t , 



dt'GtG* 



p t',Po 

Pt, [%t',P ] 

p t' , Po 
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which can be easily calculated due to the presence of 
delta-functions in Eqs. ([55]) (|57|) . resulting in the follow- 
ing expression: 



X 2 (*) 



= A t e 2Tt [x t , [x u p ] 



B, 



Pi 



+C t e 



x t , 



t,Po 



p t,Po 
Pu [xt,p ] 



(40) 



Continuing these steps, we see that the structure of the 
term \ 2 it) is repeated each time, resulting finally in re- 
placing the initial operator p by the time dependent op- 
erator p t . Thus we obtain the following closed equation 
governing the time evolution of the statistical operator 
averaged over stochastic forces (see Appendix A for de- 
tails of the derivation): 



1 

ih 



Hef(t),f> t 



j^(A t e 2Tt [x,[x,p t ] 



-B t [P,[P,p t }}+2C t e^[x 1 [P,p t }}) ■ (41) 



dp t 
dt 



We would like to emphasize that no truncations of higher 
order terms were done in deriving Eq. (|4ip . so this 
equation holds for any TD force strengths (coefficients 
A t ,B t ,C t ). A possible additional term proportional to 
[P, [x,p t ]] in the RHS of Eq. (jlTj) is redundant because 
of the identity [x, [P,p t ]] = [P, [x,p t ]]. Had we assumed 
one force F t only, the coefficients B t and Ct would be 
zero. Equation (|4T|) is structurally analogous to that ob- 
tained by Hu, Paz and Zhang [78J, although their TD 
coefficients were derived assuming a reservoir made of 
HO's, while ours are purely phenomenological. Below we 
show that the associated Wigner functions coincide, and 
as such Eq. (T4"l"j) contains the non-Markovian effects (in 
the sense of Refs. [76|, [78|, [8l|, [Hj], i.e., time-dependent 
diffusion and drift coefficients). In our derivation, limits 
as high or low-temperature, strong or weak-coupling have 
no room for discussion since all the effects of the environ- 
ment on the system depend on the adopted values for the 
four TD parameters. In particular, for coefficients A, B 
and C being time-independent one retrieves the Marko- 
vian limit. 

Eq. (T4"Tj) cannot be immediately identified with Eq. 
([6]) , due to different meanings of the variables p (physical 
momentum) in © and P (canonical momentum) in (|4ip , 
besides the presence of the factors exp(r t ) and exp(2r t ). 
However, it is easy to show that these equations are 
physically equivalent, because they give the same Wigner 
function for the mapped physical coordinates. 

The equations of motion of first and second moments 
for the canonical variables are (the average values are 
defined as (i)= Tr (Ap t )) 



d (x) /dt = e- L * <P) /to, d(P) /dt = (-dV/dx) 



d (x 2 )/dt = e~ Tt ({£, P}) /to + 2B t , 
d{P 2 )/dt = -({P,dV/dx})e Tt + 2A t e 2F \ 



4<{i,p}H-<pV 

dt to 



4C t e r < 



and one notices that these equations contain the TD 
exponential factors exp(±I\). However, passing to the 
physical momentum ([25]) , we get rid of these factors, 



d (x) /dt = (p)/m, 



d(p)/dt = -(dV/dx)-t t (p) 



d(x 2 )/dt = to" 1 {{x,p}) + 2B t , 



(42) 



(43) 



(44) 



d(f)/dt = -({p,dV/dx}) - 2t t (f) + 2A U (45) 



j f ({x,p}) - ^(p 2 )-2(x^)-t t ({x,p})+4C t . (46) 



Comparing Eqs. dHJ) and (TJ3]) with (gSJ and gSJ), we 
see that they coincide if V(x) = muf^x 2 /2, T t = 2Xt, and 
p = A. Then comparing equations for the covariances of 
the physical momentum and coordinate, one can verify 
that they satisfy the matrix equation (|14p if the diffusion 
coefficients are identified as 



At = D x (t) , B t = D p (t) , C t = D z (t) . 



(47) 



In what follows we will show the equivalence between 
the Wigner function W p (Q,t) in the canonical phase 
space, Q = (x, P) with the Wigner function in the phys- 
ical variables phase space W(q, t). The W P (Q, t) is gov- 
erned by the Fokker-Planck equation pip with time- 
dependent drift and diffusion matrices, 







,-At 



,2„Xt 



/to 







XI. 



B t C t e 
C t e xt A t e 2Xt 



The solution of Eq. (Ti~Tj) for the function VK(Q, t) can be 
written as 

W p (Q,t) = jG p {Q,Q\t)W p (Q'fl)dQ', 
where the propagator is given by the formula (40l . [52l | 
G p (Q,Q',t) = (2vr v /detN(t)" 



x exp 



i(Q - R p (t)Q')N- 1 (t)(Q - R p (t)Q') 



(48) 



Matrix N(t) satisfies Eq. pT|) (with matrices A p and 
D p ) and the initial condition N(0) = 0, whereas matrix 
R p (t) satisfies the equation dR p /dt = A p (t)R p and the 
initial condition R p (0) = 1 (unity matrix). 
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Differential equations for three different matrix ele- 
ments of the symmetrical matrix N (t) have the form 

dNn/dt = 2e- xt N 12 /m + 2B t , 

dN 12 /dt = e- xt N 22 /m-mujle xt N 11 + 2C t e xt , 

dN 22 /dt = -2mwle xt N X 2 + 2A t e 2Xt . 

Doing the substitution (|4"7]) and the changes 

N 22 = e 2Xt M 22 , N 12 = e xt M 12 , N n = M n 

we obtain the equations for elements of matrix M (t) with 
time-independent drift and diffusion matrices, given by 
Eqs. (|12p and (fTl| with A = \i. Then one can verify 
that the propagators in the canonical and physical phase 
spaces are related by a simple formula 



G P (Q,Q',t) = e- A *G(q,q',i) 



(49) 



Therefore, as the initial Wigner function is the same 
in both coordinate systems, W p (Q,0) = VF(q, 0), the 
Wigner function at any time becomes 

W p (Q,t) = J G p (Q;Q',t)W p (Q',0)dQ' 

e~ xt G(x,p; x',p',t)W(x',p';0)e Xt dx'dp' 



G(x,p; x ,p , t)W(x ,p ; 0) dx dp 

= W(q,t). (50) 

Thus the master equation (|41[) is completely equivalent to 
the translationally invariant (A = fj,) form of the master 
equation for quantum Brownian motion ^ with opera- 
tors © and ©. 

IV. SUMMARY AND CONCLUSIONS 

We have analyzed the phenomenological approach to 
build a master equation for describing the irreversible and 
dissipative dynamical evolution of the state of a quantum 
system S, under the influence of an unspecified environ- 
ment. In contradistinction to the microscopic approach 
that models the environment as a reservoir 1Z made of an 
infinite number of degrees of freedom (for example, har- 
monic oscillators), the phenomenological approach makes 
use of dissipative functions and stochastic forces. We 
showed that the Newton equation of motion for S does 
not change by introducing two such forces instead of 
only one, however when we do the quantization of the 
system, both forces become quite important, contribut- 
ing on equal footing. We derived from first principles, 
with a Lagrangian containing one dissipative function 
and two stochastic forces, the master equation describ- 
ing the quantum Brownian motion (of a harmonic oscilla- 
tor) with translationally invariant damping and the most 
general bilinear (with respect to the coordinate and mo- 
mentum operators) relaxation superoperator, which can 



be reduced to an equivalent differential equation of the 
Fokker-Planck type. However, the TD phenomenologi- 
cal parameters entering the forces cannot be determined 
within the framework of the phenomenological approach; 
they should be fixed either from experimental data that 
reproduce relevant physical properties of S or from some 
other considerations, such as the requirement of satisfy- 
ing the positivity constraint (|18[) or by fitting properties 
at thermal equilibrium. In this direction we verified that 
the master equations derived in [7(1 [73, [zH , containing 
non-Markovian effects (present in their TD coefficients) 
are accounted in our derivation where the environment is 
simulated by two effective forces and the dissipative func- 
tion, instead of assuming an interaction with an infinite 
set of HO's. 

We would like to emphasize that by averaging over 
the stochastic forces, we did not disregard any term (see 
Appendix for details). In this sense, the phenomeno- 
logical derivation of the master equation is as exact as 
other approaches [78j]. Of course, this happened due to 
the choice of correlation functions in the form of delta 
functions, albeit multiplied by time-dependent strength 
factors. The pertinent question, what could happen in 
the most general non-Markovian case, when correlation 
functions, such as A{t\,t 2 ), are arbitrary functions of the 
time difference t\ — t 2 (colored noise), requires a sepa- 
rate study. Certainly, the phenomenological approach 
used in this paper has limitations, because it is based on 
some effective Lagrangian. Therefore, although it works 
well for one-dimensional systems (or isotropical multidi- 
mensional ones), it will fail for generic multidimensional 
systems with several independent damping coefficients, 
because no effective Lagrangian can be found for such 
systems [III [86| . This explains also why only a subfam- 
ily of master equations ©, restricted by the condition 
A = /i (translationally invariant damping), can be ob- 
tained within the framework of the scheme used in this 
paper: there are no effective Lagrangians for A / fi [i.e., 
for two "friction forces" in the classical equations of mo- 
tion and ©]. 
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APPENDIX A: DERIVATION OF THE 
PHENOMENOLOGICAL MASTER EQUATION 

Iterating Eq. (|3"!2)) once more we get the formal solution 



Pt = Po + Jfr I 



W(h),p 



o 

2 ft r-ti 

dt\ I dt 2 
(i Jo 



(Al) 



W(h), W(t 2 ),pu 
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and deriving it with respect to time we get Eq. (|33[) . Us- 
ing the RHS of (|A1[) for p t , , inserting it in the RHS of Eq. 
(|33|) and keeping even terms in W (because only the even 
terms will survive after averaging over the ensemble) , Eq. 
(I33p may be written as an infinite series 



1 



dt' 



W(t), W(t'),p + 



1 



(my 



dt (tny Jo 

x / dt! f 1 dt\ \w(t ± ), Iw^rpo 



or 



dp t 
dt 



which is exact. Note that on the RHS of Eq. (|A3[) the 
density operator is time-independent whereas it is TD in 
Eq. llAll) . 

To illustrate the calculations, let us consider an exam- 
ple of quartic terms of the form 

— !-j f dt' f dt" f dt'" exp (T t + IV + IV + IV) 
(in) Jo Jo Jo 

xF t F t ,F t „F t „, [x t , [x t >, [i t //,[5 t ///,p ]]]] . 
They give the following terms in the master equation: 



(ih) Jo 
1 



dt' W(t), W(t'),p 
dt' W(t) 



dp t 



(ih) Jo 
x I dti I dt[ 



k=i w 



dt 2 e ri+r ^A t 6(t-t 2 )[xt,[x t2 ,p ]\ 



o 







dt 



dt' h 



k I UL k 




W(h), W(t[),... W(t k ), W(t' k ), Po 



(A2) 



with t' = t' . Averaging the quadratic term in W (■) of 
(|A2[) , we get the four terms of Eq. ([TO)) . Averaging the 
quartic terms in W (•), 



/ dt' [ dh [ dt[ 

(ih) Jo Jo Jo 



W(t), W(t>), W(h), W(t[),p 



we obtain sixteen terms, which together with the four 
terms of (|40|) give the expression 



dt (ih) jo 

o2 pt pti pts 

+ 4 / dt 2 / dt 3 / dUe Vt+Tt i +Tt i +Tt * 

(ih) Jo Jo Jo 

x [xt, [x t2 ,[xt 3 ,[x t4 ,p Q }}}} 

x {A t S (t - t 2 ) A t3 5 (t 3 - t 4 ) 

+A t S(t-t 3 )A t2 6(t 2 -t i ) 

+A t2 5 (t 2 - t 3 ) A t S (t - u)} + ■ ■ ■ . 

Due to the time ordering t > t 2 > t 3 > t 4 , only the prod- 
ucts of delta- functions S (t — t 2 ) S (t 3 — t 4 ) contribute to 
the integral, and in general, 



F tl F t2 ...F t2n =2 n l[A t2i _ 1 5(t2 i -i-t2i). (A5) 

i=l 

Thus the terms proportional to the coefficient A t can be 
combined as follows, 



{A t e 2r < [x (t),[x (t),p t ,} 



(ih) 

+C t te Vt (jx (t), [P(t),~Pt> 
+B t [P(t),[P(t),~Pt]]}, 



dp t 



A 



* c 2r t r~ 



dt - {lh f e [XtA ^ pM + \my\2 / 

x / dt 3 A t . 3 e 2r ^ [x t , [x t , [x t3 , [xt 3 ,Po]}}} + 
Jo 



2 2 A t (1 



P(t),[x (t),p t > 



where 
Pt' = Po 



A 



I r 

—2 / dt" {A t „e 2T t" \x t »,[*v',fa] 
[ih) Jo 



(ih) 

if one notices that 



l —e 2Tt [x t , [xup t ] 



(A6) 



+CVe I > ( 



X t " 



Pt" , Po 



-B t " 



Pt", 



Pt»,Po]]} 



Pf, fit", Po\ ) 



(A3) 



Pt = Po + 



1 



(ih) Jo 



dt 3 A t3 e '3 [x t3 ,[xt 3 ,p ]] + 



Finally, after averaging and collecting all terms in Eq. 
(|A2[) one obtains the master equation in the interaction 
picture 



dp t 
dt 



1 



2 {A t e 2Ft [x u [xt,Pt]]+B t 



(ih) 
+C t e r < ( 



x t , 



Pt,Pt 



Pt, 
Pt, [x t ,Pt] 



Pt,p t 

} , (A4) 



The factor 1/2 in the first line of (|A6j) occurs due to 
the integration of the Dirac delta function over a semi- 
infinite interval. We emphasize that Eq. (|A6|) is exact 
under the assumption (|35p , because no terms were disre- 
garded. 

Returning to the original operator p t with the aid of 
the transformation (|2T)|) , we arrive at Eq. fld]) . 
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